Stabilization of a Bose-Einstein droplet by hyperfine Rabi oscillations 
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A self-trapped Bose-Einstein condensate is shown to be stabilized in two-dimensional free space 
by Rabi oscillations between two hyperfine states which make an effective interatomic interaction 
oscillate in time. The stabilization mechanism is elucidated by using a two-component Gross- 
Pitaevskii equation combined with a variational analysis. The parameter regime of stability is 
investigated. 
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I. INTRODUCTION 

A droplet of water is a self-trapped object, in which 
attractive and repulsive interactions between the water 
molecules are balanced. Such a self-trapped system is 
hard to achieve in the usual gaseous phase, since ther- 
mal expansion of the gas must be prevented by a strong 
attractive interaction between particles, which, however, 
would lead to collapse of the gas into a denser phase. 
A gaseous Bose-Einstein condensate (BEG) is the most 
promising candidate for creating a self-trapped gas, or a 
"BEG droplet," because of unprecedented controllability 
of the parameters of the system. 

In one dimension (ID), a BEG droplet, or a matter- 
wave bright soliton, is stable with an attractive inter- 
action between atoms counterbalancing zero-point quan- 
tum pressure. This novel state of a gas has been exper- 
imentally achieved by the ENS group [l| and the Rice 
group [2] by using a ''Li condensate in quasi-lD traps. 
In 2D and higher dimensions, however, there is neither 
a stable nor metastable state of a BEG droplet, if the 
interaction is short-range and constant in time 3] . 

Dynamic stabilization is a possible way to obtain a 
BEG droplet in 2D and 3D. It has been shown that if the 
interaction is made to oscillate in time, a BEG droplet 
can be dynamically stabilized in 2D 4, ^] by a mechanism 
similar to that for the stabilization of an inverted pendu- 
lum . This stabilization mechanism for a BEG droplet 
has been studied by several authors 0, H, H, [13, [0, 13 ■ 
We have shown that the same mechanism can also sta- 
bilize a BEG droplet in 3D in the presence of dissipa- 
tion [l^. Another dynamic scheme to stabilize a BEG 
droplet_in 3D is to use feedback control of the interac- 
tion 
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In respect of the stabilization method proposed in 
Refs. [1, we must oscillate the interaction between 
repulsive and attractive at a frequency much higher than 
the characteristic frequencies of the system. If we use the 
magnetic Feshbach resonance [is'] to oscillate the interac- 
tion, we should oscillate the strength of the applied mag- 
netic field at such a high frequency. We propose in the 
present study a method to achieve this dynamic stabiliza- 



tion without using any time-dependent Feshbach control 
of the interaction. Our idea is to use two hyperfine states 
with different scattering lengths. We will show that the 
Rabi oscillation between these hyperfine states causes an 
effective oscillation of the scattering length, thereby sta- 
bilizing a BEG droplet through the mechanism proposed 
in Refs. [1, . This method simply irradiates the system 
with a constant electromagnetic wave. 

This paper is organized as follows. Section |TT] is de- 
voted to describing a system of a two-component BEG 
driven by an external field. Section [llll performs numer- 
ical integration of the two-component Gross-Pitaevskii 
equation to show that the BEG droplet is indeed dynam- 
ically stabilized. The stability diagram for the scattering 
lengths is obtained. Section llVl derives an effective single- 
component equation, and combines it with a variational 
method to elucidate the stabilization mechanism for a 
BEG droplet by hyperfine Rabi oscillations. Section Ivl 
provides conclusions. 



II. FORMULATION OF THE PROBLEM 

We consider a situation in which bosonic atoms with 
hyperfine degrees of freedom are irradiated by a uniform 
electromagnetic wave that is resonant with the energy 
difference between two hyperfine states of the atoms. The 
atoms are assumed to be tightly confined in the z direc- 
tion by harmonic potential raLo^z^ jl that is independent 
of the hyperfine state. If hio^ is much larger than the 
other characteristic energy scales, the wave function in 
the z direction is frozen in the ground state of the har- 
monic potential, and then the system can be reduced to 
2D. We will consider the 2D system in the subsequent 
analysis. 

The Hamiltonian for the system consists of three parts: 



int ■ 



(1) 



The non-interacting part of the Hamiltonian is given 
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by 
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Mr), (2) 



where annihilates an atom in hyperfine state i, m 

is the atomic mass, and Vi{r) is an external trapping 
potential for hyperfine state i. The effect of the resonant 
field on atomic state is described by 



Hf = h dr 



(3) 



where jfij is the Rabi frequency and arg(ri) is the phase 
of the field. The interaction between ultracold atoms is 
short-range and described by the contact Hamiltonian, 
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dr ^ ay-(/>|(r)^](r)t/>j(r)t/>,(r), 

(4) 

where a^j is the s-wave scattering length between atoms 
in hyperfine states i and j. The factor [mujz/ {2'kK)Y/'^ in 
Eq. ^ comes from the integration with respect to z. 

In the mean-field approximation, the system at zero 
temperature is described by the 2D Gross-Pitaevskii 
(GP) equation, 

ot \ 2m I 



+5l2|V'2pV'l, 



(5a) 



ih^^ - ( -:^V^+y2 1^2+ 522 1^-2 1'^2 



dt 



2m 

-5l2|V'lPV'2, 



(5b) 



where V'i('"i i) is the macroscopic wave function satisfying 
the normalization condition, 



/ dr^|V.(r,t)|2 = Ar, 

i=i 



with N being the number of atoms, and 
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(6) 
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TO V 2Trh 
representing the interaction coefficients. 



III. STABILIZATION OF A TWO-COMPONENT 
BEC DROPLET 

The interaction parameter of the system oscillates be- 
tween gii and 1722 by the Rabi oscillation between two 
hyperfine states. If one of them is positive and the other 
is negative, the stabilization mechanism in Refs. [1, [H] is 
expected to be applicable to the present system. This is 
the main idea of the present paper. 




FIG. 1: (Color) (a) Time evolution of the normalized peak 
density di+d2 and that of the fraction of atoms within r = 100 
fim for an = 0.75 nm, 022 = —1.5 nm, ai2 = 0, and f2 = 27r x 
500 Hz. The initial state is the ground state in the presence 
of a radial trapping potential with u!± — 2-k x 25 Hz, and the 
trapping potential is switched off at t = 0. (b) Time scale 
of (a) is magnified. Small- amplitude oscillations in (b) are 
invisible in (a) because of the time resolution, (c) Snapshots 
of the normalized density profiles (IV'iP + |^/'2 |^)?i/(?Ti<^±A'") 
at t = 0, 430 ms, and 500 ms. The size of the images is 30 
fim X 30 ^m. 



We assume that all the atoms are initially in hyper- 
fine state 1, and confined in radial trapping potential 
Vi ~ miJ^r^ 12. The initial wave function V'l is then the 
ground state of GP equation ((5a| with ^^2 = and SI = 0. 
At < = 0, the radial trapping potential is switched off and 
the radiation field is switched on, i.e., V\ = V2 = and 
17 7^ for t > 0. 

We first investigate dynamics of the system for specific 
values of scattering lengths: a\\ = 0.75 nm, 022 = —1.5 
nm, and an = 0. The radial trapping frequency for 
preparing the initial state is w_l = 27r x 25 Hz, and the 
trapping frequency for confinement in the z direction is 
tJz = 27r X 5 kHz. The ratio lOz/lo^ — 200 is sufficient for 
the 2D approximation to be valid. The number of atoms 
is iV = 2500, and the Rabi frequency is = 27r x 500 
Hz. The ratio w^/fi = 10 assures that no axial modes 
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are excited by the Rabi oscillation. For this condition, 
the attractive interaction by 022 dominates the kinetic 
pressure, and the system collapses even in the absence of 
the trapping potential if the atomic state is fixed to the 
hyperfine state 2. 

Figure [T] (a) shows the time evolution of normalized 
peak density di + d2 and that of the fraction of atoms 
around the center, Nioofj.m{t)/N, where 



d^{t) = 



muj±N 



|^,(r = 0,t)p (i = l,2), (8) 



(9) 



We can clearly see that collapse is prevented and the 
BEC droplet is dynamically stabilized. The Rabi oscilla- 
tion between the two hyperfine states induces oscillations 
of the total density at the same frequency, with the am- 
plitude being larger for higher density (compare the left 
and right panels in Fig. [1] (b)). The long-period oscilla- 
tion at the frequency of about 8 Hz persists for a long 
time (solid curve in Fig. [T] (a)). The initial decrease in 
the dashed curve in Fig.[T](a) shows that a small fraction 
of atoms are lost from the central region mainly during 
the initial formation of the droplet. These atoms have 
escaped from the droplet for lack of the radial trapping 
potential. 

We next consider the stability of the BEC droplet for 
various values of scattering lengths. Figure [2] (a) shows 
stability diagram with respect to an and 012 for an = 
— 1.5 nm, and Fig. [2] (b) with respect to a22 and ai2 
for an = 0.5 nm. In the white region, the lifetime of 
the BEC droplet is at least 2 seconds. We find that a 
band of stable region lies between regions of collapse and 
expansion. These two types of instability are divided by 
the solid lines in Fig. [2] The stable region disappears 
for ai2 ^ —0.2 nm. As we can see in Fig. [5] (a), the 
value of an in the stable region become negative as ai2 
increases. However, this would not hinder an initial state 
of the condensate from being prepared. One can prepare 
the condensate in a positive an and switch an to an 
appropriate negative value at t = 0. 

The stability region is not sensitive to the Rabi fre- 
quency fl. To see this, we normalize the time, length, and 
wave function by \n\-^, [h/{m\n\)]^/^, N[m\n\/h]^/^, re- 
spectively. CP equation (|5a|) then becomes 



2 

-|-47rA^ 



2\n\ 



27rh 



(anlV'iPV'i + ai2|V'2pV'i) 



(10) 



We note that fl is included only in the term of the radial 
trapping potential used for preparing the initial state, 
and therefore the change in Q only affects the initial wave 
function in the normahzed form of CP equation pU)) . We 
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FIG. 2: Stability diagrams for a BEC droplet with respect 
to (a) ai2 and an for 022 = —1.5 nm, and (b) ai2 and 022 for 
an = 0.5 nm. The gray scale shows the lifetime of the BEC 
droplet. The white region has a lifetime of at least 2 s. The 
solid lines divide the regions of collapse and expansion. The 
dashed lines are given by Eq. 



performed numerical simulations also for Q = 2ttx1 kHz, 
and confirmed that the stability region is almost the same 
as in Fig. 



IV. ANALYTIC RESULTS 

In this section, we analytically study the stabilization 
of a BEC droplet shown in Sec. IIIII 

We first transform wave functions ijji and "02 into 

ifi = ■il;icos\n\t + ie"''il;2sm\n\t, (11a) 
*2 = V'2Cos|17|i + ie"*"*V'isin|rj|t, (lib) 

where = e''^|rj|. Substituting Eq. ^ into CP equa- 
tion (O, we find that the terms proportional to jfij van- 
ish, and we obtain 
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ot 2m 8 

+ (311 - 522) |^^^|*iP*i + [e-*^vl>;.l/2 _ (2|vi.,|2 + |*^|2) | 

8 

+i sin4|ri|i [e-^^*^*? + e'^ (-2|*i|2 + l^-sP) ^2] }, (12) 



where we set Vi — V2 — 0. The expression of ifid'i>2/dt 
is obtained by replacements ^'i ^ 5*2, 511 522, and 



^ e 



in Eq. ((T 



The initial state of V'2 is assumed to be ^"2 = and 
then ^1/2 = at t = 0. Setting ^'2 = in the expression 
of ihd^2/dt gives 



ih 



dt 



— [2{gn- 922)sin2\n\t 
<9ii + 922 - 2512) sin4|r!|t] |*i|2*i,(13) 



which indicates that ^'2 is of the order of l/|ri|. We there- 
fore approximate that ^'2 is always zero. Equation \\2\ 
is then approximately reduced to the single-component 
GP equation, 



ih 



^ ^2 T 



G(t)|*i|'*i, 



(14) 



where 



with 



G(t) =G'o + Gicos2|r2|i-f G2Cos4|r2|i, (15) 



Go = 
Gi = 
G2 = 



3(ffii +522) + 2512 



511 - 522 



.gii + 522 - 2gi2 



(16) 
(17) 
(18) 



To check if Eq. captures the stabilization mecha- 
nism for a BEG droplet, we perform a variational analysis 
by using a Gaussian variational wave function ^16.] , 



*i 



exp 



2i?2(t) ' " 2hR{t) 



+ i 



(19) 



where R(t) is the variational parameter that character- 
izes the size of the droplet, and the second term in the 
exponent describes mass current which is required to sat- 
isfy the equation of continuity. The action that derives 
Eq. (fT4|) has the form. 



K 



drdt 



(20) 



Substituting Eq. (HH) into Eq. ^ and taking SK/SR = 
0, we obtain the equation of motion for R{t), 



Rit) = 



NG{t) 



27rm 



R'it)' 



(21) 



Noting that the dynamics comprises rapid and slow oscil- 
lations as shown in Figs. [1] (a) and[l](b), we separate R{t) 
into a slowly varying part Ro(t) and rapidly oscillating 
part e{t) as 



R{t) = Ro{t) + e{t). 



(22) 



The rapidly oscillating part of Eq. (PT|) is approximated 
to be 



N 



2TTmR^{t) 



(GiCos2|17|t-)-G2Cos4|r2|t), (23) 



where we neglect the terms of order Since R^ 

can be regarded as a constant in the fast time scale in 
Eq. ((231), we have 



sit) - 



N 



327rmi?3(<)|f7|2 



(4Gi cos2|r2|i + G2 cos4|17|i) , 



(24) 

Substituting Eq. (|24|) into the slowly varying part of 

Eq. m, 



Ro{t) 



^ NGo] 1 
m2 2TTm J Rl{t) 
3N 



2TTmR^{t) 
+0(e^), 



e{t) (Gi cos2|l]|t + G2 cos4|r2|t) 

(25) 



where indicates time average for the rapid oscilla- 
tion, we finally obtain an effective equation of motion for 
slowly varying parameter Ro{t) as 



Ro{t) ^ 



NGq \ 1 3iV^ (4G^ + Gj) 
27rm J i?3(t) ^ 128TTm^\n\^Rlit)' 

(26) 

The first term on the right-hand side of Eq. orig- 
inates from the kinetic energy and the constant part of 
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the interaction energy. For stationary Rq to exist, the 
coefficient of the first term must be negative, i.e.. 



mNGo 
2nh^ 



< -1, 



(27) 



which gives the condition for stabihty against expansion. 
Using Eqs. (fTHll and ((27|) . the boundary between the sta- 
ble and unstable regions is given by 



3(aii + 022) + 2ai2 



4 

'n 



2Trh 



mojz 



(28) 



This boundary is shown in Fig. [2] by the dashed lines, 
which show close agreement with the numerically ob- 
tained boundary. 

We note that the second term on the right-hand side 
of Eq. ([26l) . which originates from the oscillating part of 
the interaction in Eq. (jl4p , prevents the system from col- 
lapsing. Thus, the effective oscillation of the interaction 
by the Rabi oscillation dynamically stabilizes the BEC 
droplet. The system is stationary {Rq — 0) with the size 
of the droplet. 



647r|f7|2 27rr 



iNGn 



1/4 



(29) 



which is ~ 6.4 /im for the parameters in Fig. [T] Compar- 
ing the density profile at t — [R = {fi/muj±_Y^'^ ~ 7.6 
/im] with those ai t = 430 and 500 ms in Fig. [T] (c), we 
can see that this value of Rq agrees well with the numer- 
ical one. The frequency of the small oscillation around 
Eq. ((29|) is obtained as 



8%/2|f7| |27rn^ -hmiVGol 
^/imN yfAG'i + G'i 



(30) 



which is ~ 27r X 12 Hz for the parameters in Fig. [U in 
qualitative agreement with the numerical result of about 
27r X 8 Hz. The difference may be due to the large am- 
plitude of the oscillation as shown in Fig[T] (a). 



V. CONCLUSIONS AND DISCUSSION 

We have studied a 2D BEC of atoms in two hyper- 
fine states, which undergo Rabi oscillation by an external 
field. When scattering lengths an, 022, and ai2 are dif- 
ferent, the Rabi oscillation leads to effective oscillation of 
the interaction. If the scattering lengths fall in the stable 
regions such as those in Fig. [21 the self-trapped BEC in 
free space — a BEC droplet — can be stabilized by the 
mechanism discussed in Ref . 0, 0] ■ 

Solving the two-component GP equation ^ numeri- 
cally, we have shown that the stabilization mechanism 
works and the BEC droplet is stabilized (Fig. [T]). The 
stabilization is shown to be possible even if the system 
starts from the ground state in a radial trapping poten- 
tial, which is switched off at i = 0. We have obtained 



the stability diagram with respect to scattering lengths 
(Fig. [5]), in which there exists the stability region between 
the regions of collapse and expansion. 

To elucidate the mechanism for the dynamic stabiliza- 
tion of a BEC droplet, we approximately reduced the 
two-component GP equation to a single component equa- 
tion ()14|) and studied it by using the Gaussian variational 
method. We found that the Rabi oscillation generates an 
effective potential that counteracts the attractive interac- 
tion, thereby creating a stable droplet state. The size of 
the droplet and its breathing-mode oscillation frequency 
are obtained by the variational analysis, which produces 
results in qualitative agreement with the numerical ones. 

Our scheme applies to any system where atoms in the 
two relevant hyperfine states are not transferred to other 
states by inelastic collisions and the coupling field is avail- 
able that makes Rabi oscillations between the two states. 
For example, the hyperfine states \F,mp) = |1,1) and 
|2,2) of *^Rb and ^^Na, and |2,2) and |3,3) of ^^Rb 
are possible candidates, if their scattering lengths can 
be tuned to fall in the stable region. Unfortunately, for 
1 1,1) and |2,2) hyperfine states of ^Li, the scattering 
lengths are an ~ 0.26 nm 0], a22 — —1.4 nm [l3|, and 
ai2 ~ — 1.4 nm [l^, and this set of scattering lengths is 
located out of the stability region. By systematic calcu- 
lations as in Fig. ^ we found that there is no stable set 
of an and a22 for this value of ai2. 

Using the above pairs of hyperfine states, inelastic de- 
cay of the upper hyperfine state occurs only through 
dipole-dipole interaction. The dipolar decay rate K is 
typically less than 10^^^ cm^/s 19] and the peak den- 
sity in Fig.[T](a) is n ~ 0.4:NmLj_i/h[mLJz/ (t^^)]^^^, which 
gives nK < 0.2 s^^. Therefore, the lifetime of the droplet 
in the presence of the dipolar loss is estimated to be more 
than 5 s. 

Another possibility is the use of two states in the 
same hyperfine manifold, e.g., 11,1) and 11,0), in which 
modest magnetic fields (~ 1-10 G) should be applied to 
suppress the spin-exchange colHsions (]!, 0) -I- Jl, 0) 
11,1) -I- 11,-1)). For example, for these states of "^Li, 
there are regions in which Eq. ([28]) becomes small and 
negative as a function of magnetic field, suggesting that 
the stabilization is possible. 

The 2D potential may not be perfectly flat in experi- 
ments, and a residual weak potential in 2D, e.g., a weak 
harmonic potential, may exist. Even in this case, we can 
distinguish self trapping from trapping by a residual po- 
tential, if the size of a droplet is made much smaller than 
the scale of inhomogeneity in the 2D potential. 

As in the case of oscillating the scattering length 0, , 
we have not been able to create a stable BEC droplet in 
3D in the absence of dissipation. One reason for difficulty 
is that collapse is stronger for higher dimensions [iS] , and 
therefore the amplitude of the rapid oscillation should 
be larger to suppress the collapse. However, this breaks 
the small- amplitude and Gaussian approximations, and 
unstable modes are dynamically induced. 

The present method is different from the one in 
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Ref. [3, 's'l , in that explicit oscillation of the interaction 
using the Feshbach resonance is not needed. The rapid 
modulation of the applied magnetic field for the Fesh- 
bach resonance may cause some experimental difficulties, 
e.g., an induced current in the apparatus and resultant 
heating, and therefore the present method is expected to 
facilitate the experimental realization of a BEC droplet. 
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